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Abstract. The Birmingham Solar-Oscillations Network (BiSON) has been 
collecting data for over 30yrs and so observations span nearly three llyr so- 
lar activity cycles. This allows us to address important questions concerning 
the solar cycle and its effect on solar oscillations, such as: how consistent is 
the acoustic behaviour from one cycle to the next? We have used the p-mode 
frequencies observed in BiSON data from one solar activity cycle (cycle 22) to 
predict the mode frequencies that were observed in the next activity cycle (cycle 
23). Some bias in the predicted frequencies was observed when short 108 d time 
series were used to make the predictions. We also found that the accuracy of 
the predictions was dependent on which activity proxy was used to make the 
predictions and on the length of the relevant time series. 



1. Introduction 

The frequencies of solar oscillations vary with solar activity ( Woodard &: Novesl 



19851 ) and the changes in frequency can, in principle, be used to infer solar-cycle 



related structural changes in the Sun. Methods have been developed to remove 
solar cycle frequency shifts (e.g. IChaplin et aT1l2004l ) to allow comparisons be- 
tween the mode frequencies obtained from non-contemporaneous data and to 
enable accurate activity-independent structural inversions to be performed. 

The Birmingham Solar-Oscillations Network (BiSON) has been recording 
data for over 30yrs, and so observations span nearly three solar cycles. We 
aim to take advantage of this to answer the following question: can the mode 
frequencies observed in cycle 22 be used to predict the mode frequencies observed 
in cycle 23? In other words, is it possible to produce a model that can be used 
to predict the mode frequencies observed in cycle 23? 

Each llyr solar cycle is different from the previous one, for example, cy- 
cle 23 has been descr ibed as magnetically simpler than recent previous cycles 
( de Toma et al. 20041 ) . Several different activity proxies can be used to mea- 



sure the level of solar activity and, as each proxy responds to different physical 
processes, the amplitude and shape of each activity cycle depends upon which 
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Figure 1. The \\ observed at different activity levels when four different 
activity proxies were used to predict the frequency shifts. 

activity proxy is considered. With this in mind we aim to determine whether 
the acoustic behaviour of the Sun is predictable from one cycle to the next. 

We begin this paper by describing how we used the mode frequencies, ob- 
served in cycle 22, to predict the mode frequencies observed in cycle 23 (Section 
2). Then, in Section 3, we compare these predicted frequencies with the actual 
frequencies observed in cycle 23. Section 4 comprises a discussion of the results. 



2. Predicting frequency shifts 

An llyr BiSON time series is available, which covers the whole of cycle 22. 
To allow solar cycle frequency shifts to be examined this data set was split 
into contiguous 108 d time series. A standard likelihood maximization metho d 
was used to fit the power spectra of these time series (e.g. Chaplin et al.|[l999l ). 



enabling the observed mode frequencies to be determined. The weighted mean 
activity level, a, observed during each 108 d time period was also calculated for 
four activity proxies, namely the 10.7cm flux (F10.7), the Hel equivalent width 
(Hel EW), the Kitt Peak Magnetic Index (KPMI) and the International Sunspot 
Number (ISN). 

A minimum activity reference set was determined by averaging the frequen- 
cies from five 108 d time series that were observed during the minimum activity 
epoch at the boundary between cycle 22 and cycle 23. We define the solar cycle 
frequency shifts as the difference between the frequencies given in the minimum 
activity reference set and the frequencies of the corresponding modes observed 
at different epochs and, consequently, different levels of activity. The size of 
a frequency shift has a well-known dependency on frequency an d mode inertia 



and th ese dependencies were removed in the manner described in lChaplin et al 



( 20041 ) . We have, therefore, determined the mean frequency shift at a given 



activity level. We now describe how these shifts were used to predict the mode 
frequencies observed in cycle 23. 
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Tabic 1. The weighted mean x"i f° r the different activity proxies 



Proxy 


Weighted mean x 2 


Distance from x 2 = 1 


F10.7 


1.18 ±0.06 


3.1(7 


Hel EW 


1.20 ±0.06 


3.5a 


KPMI 


1.30 ±0.06 


5.2a 


ISN 


1.18 ±0.06 


3.1(7 



For low-degree p modes, the relationship between activity and mean fre- 
quency shift is approximately linear. A linear least squares fit was performed to 
determine the coefficients of the relationship between the mean frequency shifts 
observed during cycle 22 and the corresponding mean activity observed by each 
of the four activity proxies. These linear relationships were then used to predict 
the size of the mean frequency shift at a given activity level, which corresponds 
to a particular time period from cycle 23. 

Let bv\ (a) be the mean frequency shift predicted by the linear relationship 
for a given activity level, a. To determine the frequency shift for a particular 
mode, 8v n i(a), with degree, I, and radial order, n, it was necessary to scale 
5v-[ (a) with mode fr equency and inertia. This was done using equation 4 of 
IChaplin et al.l (|2004l ). The resulting frequency shift, 5v n ^{a), was added to the 
corresponding mode frequency stored in the minimum activity reference set to 
give a predicted frequency, and its associated errors, for a given epoch. The 
predicted frequencies were compared to the frequencies observed in BiSON data 
during the same time span and we now go on to describe the results of this 
comparison. 



3. Analysis of results 

To determine how accurate the prediction method is we have compared the pre- 
dicted frequency of a mode at an activity level, a, with the frequency observed 
in BiSON data at the same activity level. The observed frequencies were de- 
termined by fitting the power spectrum of the corre sponding 108 d time series 
using the standard likelihood maximization method ( Chaplin et al.|[l999l ). 



The reduced x between the predicted and the observed frequencies at a 
particular activity level was determined. Let x\ De the reduced x 2 f° r an modes 
for which accurate fits to the data were obtained in a 108 d time series, observed 
during cycle 23. Errors on the reduced x 2 were taken to be ^/2/{N - 1), where 
./V is the number of data points used to calculate the reduced x 2 ■ Figure [Q shows 
the variation of x\ with activity, where each panel shows the results predicted 
by a different activity proxy. The vast majority of the observed x 2 are within 
3cr of unity, which is the expected value if the model is accurate. However, for 
all activity proxies, the observed x\ are predominately greater than unity, which 
implies there is some bias in the predicted results. 

To determine whether there is any correlation between the level of activity 
and the accuracy of the predicted frequencies the Pearson's correlation coefficient 
between activity level and x\ was calculated (see the values above the panels in 
Figure [1]). None of the correlation coefficients are significant at a 95% confidence 
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Hel EW 




Figure 2. The \\ i observed at different radial orders, n, when four different 
activity proxies were used to predict the frequencies. 



level. However, the correlation coefficient for the 10.7 cm flux result is smaller 
than for the other activity proxies. The results imply that the accuracy of 
the predictions made using the Hel EW, the ISN and the KPMI are marginally 
dependent on activity. We propose that this is true for the following reason. The 
Hel EW is mainly sensitive to the weak magnetic field and the ISN and KPMI 
are most sensitive to the strong magnetic field, while the 10.7 cm flux shows 
good sensitivity to both the strong and the weak magnetic flux. Therefore, we 
believe that to make accurate predictions at all activity levels it is necessary for 
the activity proxy to be responsive to both the strong and weak flux. 

As a measure of the overall accuracy of the predictions the weighted mean 
Xa was determined, and the results are shown in Table [TJ All of the weighted 
mean x\ are more than 3a from unity, indicating there is some bias in the 
predicted frequencies. The weighted mean x\ 1S noticeably larger for the KPMI 
than when the other proxies are used, which highlights the inadequacies of this 
proxy as a predictive tool. 

We also investigated whether the accuracy of the predicted mode frequencies 
was dependent on frequency (or n) . Let x\ i be the reduced x 2 f° r individual 
n and I across all activity levels. The results are shown in Figure [2j All of the 
proxies enable the frequencies to be predicted to a reasonable degree of accuracy 

108d 216d 432d 
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Figure 3. The variation of \a with activity for time series of different 
lengths, when 108 d time series were used as the reference set. 
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Figure 4. The variation of x„ i with n for time series of different lengths, 
when 108 d time scries were used as the reference set. 



as Xn i i s within 3<r of unity for the vast majority of n and I. Each panel of Figure 
[2] contains an outlier, which is observed because the frequency determined by 
the fitting procedure is inaccurate since at certain times this mode contains little 
power. 

Although not shown here, the Pearson's correlation coefficients between 
i and n were found to be approximately zero implying the accuracy of the 
predictions is not dependent on n. However, it should be noted that this is only 
true for the range of n examined here. Once we move outside this frequency 
range (2200 fjHz < v < 3500 /iHz) the quality of both the fitted frequencies and 
the predictions deteriorates significantly. 

3.1. Time series of different lengths 

Next, we sought to determine whether the data observed in 108 d time series 
in cycle 22 could be used to predict the acoustic behaviour observed in longer 
(216 d, 432 d) data sets from cycle 23. Figure shows the variation of x\ with 
activity when predictions were made for different length time series. The deter- 
mined Xa increases with the length of the time series, implying that the predic- 
tions for longer time series are less accurate. Additionally, the weighted mean 
Xa> which are quoted in Table [2] (when the 108 d time series are used to make 
the predictions), are found to be larger for the longer time series, indicating the 
108 d time series can not be used to predict the frequencies observed in longer 
time series. 

Figure U] shows the variation of xf n with radial order, n, for different length 
time series. The determined xf n increase as n decreases, especially in the longer 
time series, indicating that the predictions for low-n modes are less accurate 
in long time series. Quoted above each panel in Figure 0] are the Pearson's 
correlation coefficients, which become increasingly negative for the long time 
series and is in fact significant at a 95% confidence level for the 432 d time series. 
This is probably because it is more difficult to observe low-frequency modes in 
108 d time series than in longer data sets and so the predicted frequencies are 
biased. 

Finally we ask whether more accurate predictions are made if the frequen- 
cies observed in longer time series are used to define the relationship between 
frequency shift and activity. For example, we have used the mode frequencies 
observed in time series of length 216 d during cycle 22 to make predictions of 
the mode frequencies observed in 216 d time series in cycle 23. The observed 
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Tabic 2. The weighted mean x\ between the observed and predicted fre- 
quency shifts for different length time series. 



Time series length Weighted mean xt Distance from xt = 1 



108 d to predict 108 d 


1.18 ±0.06 


3.1a 


108 d to predict 216 d 


1.24 ±0.07 


3.5cr 


216 d to predict 216 d 


1.06 ±0.07 


0.9cr 


108 d to predict 432 d 


1.51 ±0.10 


5.3cr 


216 d to predict 432 d 


1.13 ±0.10 


1.4a 


432 d to predict 432 d 


1.07 ±0.10 


0.8cr 



weighted mean x\i quoted in Table [21 indicate that the most accurate predic- 
tions are made when the time series used to make the predictions and the time 
series for which the predictions are required are both of the same length. 

4. Discussion 

We have shown that, when 108 d time series are used, there may be some evidence 
for a bias in the predictions, as the weighted mean Xa are noticeably greater 
than unity, indicating that it may be problematic to use the 108 d time series to 
predict the mode frequencies. The accuracy of the predictions made using the 
10.7 cm flux show least dependence on activity, probably because the 10.7 cm 
flux is sensitive to both the strong and the weak magnetic field. The predictions 
made using the KPMI were not as accurate as those made using the other models 
and so this proxy should not be used as a predictive tool. 

The accuracy of the predictions decreases when 108 d time series are used 
to predict mode frequencies for longer time series and so 108 d time series can 
not be used to predict the acoustic behaviour observed in longer time series. 
The predicted frequencies were more accurate when longer data sets were used 
to make the predictions and when the time series from cycle 22, which were used 
to make the predictions, and the cycle 23 time series were of the same length. 
Note that this analysis has been performed using the mode frequencies observed 
in 2 solar cycles only but would ideally use the data from more solar cycles. 
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